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ABSTRACT 

Context. The gamma-ray binary LS I +61°303 shows multiple periodicities. The timing analysis of 6.7 yr of GBI radio data and of 6 
yr of Fermi- LAT GeV gamma-ray data both have found two close periodicities /Lgbi = 26.49 ± 0.07 d, P 2 ,gbi = 26.92 ± 0.07 d and 
P Uy = 26.48 ± 0.08 d, P Xy = 26.99 ± 0.08 d. 

Aims. The system LS I +61°303 is the object of several continuous monitoring programs at low and high energies. The frequency 
difference between and v 2 of only 0.0006 d _1 requires long-term monitoring because the frequency resolution in timing analysis is 
related to the inverse of the overall time interval. The Owens Valley Radio Observatory (OVRO) 40 m telescope has been monitoring 
the source at 15 GHz for five years and overlaps with Fermi-LNT monitoring. The aim of this work is to establish whether the two 
frequencies are also resolved in the OVRO monitoring. 

Methods. We analysed OVRO data with the Lomb-Scargle method. We also updated the timing analysis of Fermi- LAT observations. 
Results. The periodograms of OVRO data confirm the two periodicities Rrovro = 26.5 ± 0.1 d and /L.ovro = 26.9 ± 0.1 d. 
Conclusions. The three indipendent measurements of P\ and P 2 with GBI, OVRO, and Fermi- LAT observations confirm that the 
periodicities are permanent features of the system LS I +61°303. The similar behaviours of the emission at high (GeV) and low (radio) 
energy when the compact object in LS I +61°303 is toward apastron suggest that the emission is caused by the same periodically (Pi) 
ejected population of electrons in a precessing (P 2 ) jet. 
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1. Introduction 

The stellar system LS I +61°303 is formed by a co mpac t object 
and a Be star in an eccentric orbit (e — 0.72+0.1 5, Casares et al. I 
20051) . Fermi-LXT observations (lAbdo et al. 1 120091) reveal a 
large outburst toward periastron (O per iastron = 0.23 ± 0.02, 
Casares et al. 2005[) along with a smaller outburst toward apas¬ 
tron (iJaron & Massi 1 120141) . Radio observations of the system 
show a large radio outbur st toward apastron exhibiting a mod¬ 
ulation of ~ 1667 days (Paredesetal.lll990t iGresorvI 120021: 
Massi & Kaufman Bernado 2009, and here Fig.l). 

Two distinct outbursts along the orbit are well understood in 
the context of accretion along an eccentric orbit (iBondi & Hovlei 
|1944|) . The expression of the accretion rate is 


M = 


4^Pwind(GM X )“ 


(1) 


rel 


where p w ; n d is the density, and v re i is the relative velocity between 
the compact object and the Be equatorial wind. As shown by sev¬ 
eral authors for the eccentric orbit of LS I +61303, two peaks re¬ 
sult: One peak corresponds to the periastron passage because of 
the highest density; the second peak occurs in the phase interval 
(toward apastron) where the reduced velocity of the compact ob- 
ject compensates for the decrease in density (jTay lorjflaLl 1992t 
(Marti & Paredesll 19951: iBosch-Ramon et allfeOObl: iRomero et al] 


12007 ). At each accretion peak, matter is assumed to be ejected 
outward in two jets perpendicular to the accretion disk plane 
as for microquasars, but for the two ejections different ener¬ 
getic losses for the particles occur. Ejected relativistic electrons 
around periastron suffer severe i nverse Compton losses bec ause 
of the proximity to the Be star dBosch-Ramon et aD 1 2006 ): the 
electrons upscatter stellar photons to higher energies, and by 
doing so, they loose their energy and become unable to gener¬ 
ate synchrotron emission in the radio band. The consequence 
is that around periastron a high-energy outburst caused by in¬ 
verse Compton process results, b ut no or negligible radio emis¬ 
sion dBosch-Ramon et akl 12006 ). This hist predicted gamma- 
ray outburst has been confirmed by Fermi -LAT observations 


(lAbdo et al, 1120091 Hadasch et al.ll2012l; lAckermann et al,lF20 1 3fc 

IJaron & Massi 1120141) . During the second accretion peak, the 
compact object is much farther away from the Be star, and in¬ 
verse Compton losses are lower: The electrons can propagate 
out of the orbital plane and generate synchrotron emission in 
the radio band and a minor gamma-ray outburst (see the pre¬ 
dicted outbursts in Fig. 2 of iBosch-Ramon et al .1 2006). Indeed, 
as quoted above, radio and gamma-ray observations confirm 
the predicted radio and minor gamma-ray outburst toward apas¬ 
tron. The spectrum of the radio outburst has been reported to be 
flat from 1.5 t o 22 G Hz by iGregorv et al.1 dl979l) . and recently, 
IZimmermannl (1201 3b observed a flat spectrum up to 33 GHz. 
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Fig. 1. Left: OVRO data vs MJD. Centre: OVRO (filled squares, black) data and GBI 8 GHz data (empty squares, red) vs ©; 0 is the fractional part 
of (f - to)/1667 with to=JD2443366.775 1Gregory! 20021) . Right OVRO (black) data and GBI 8 GHz data (red) vs ®, where <T> is the orbital phase, 
i.e., the fractional part of (t - t 0 )/P{, with P j = 26.496 days dGregorvll2002l) . The largest radio outbursts occur around ® = 0.6. At ® = 0.23, i.e., 
at periastron, only a low level of emission is present. 


These radio observations corroborate the microquasar hypothe¬ 
sis for LS I +61°303 because a flat radio spectrum is indeed as- 
sociated with the conical radio je t of microquasars (iFender et al.1 
2000: lFenderi200n lKaiseHl2006l) . 


In 1993, the radio-emitting source was resolved with high- 
resolution radio observations showing a structure of milliarc- 
second (mas) size corresponding to a few AU at the dis¬ 
tance of 2.0 kpc (lMassi et al. 1119931) . Successive observations 
revealed that the radio morphology not only changes posi¬ 
tion angle, but it is even sometimes o n e-sided and at other 
times two-sided ( Peracaula et a . 1 119981: iParedes et akl Il998t 


iTavlor et al. 1 12000 : Massi et al. I 2001 : lMassi et al J 120041) . This 
suggested that LS I +61°303 might be a precessin g microquasar 
(Kaufman Bernado et al.ll2002t lMassi et al.ll2004l) . A precession 
of the jet leads to a variation in the angle between the jet and 
the line of sight and therefore to variable Doppler boosting. 
The result is a continuous variation in both the position an¬ 
gle of the radio-emitting structure and its flux density (lMassi I 
12007 ). Concerning the precession time-scale, rather fast posi¬ 
tion angle variations of almost 60°/day measured by MERLIN 
observation s (lMassi et al.1 2004 ) were confirmed by VLBA ob¬ 
servations (iDhawan et alJ 2006), which measured a rotation of 
roughly 5° - 7° in 2.5 hrs (i.e., ~ 60°/day). Radiojistrometry re¬ 
sulted in a precessing period of 27-28 days (lMassi et alj|2012l) . 

A timing analysis of the GBI radio data of LS I +61°303 has 
revealed two rather close frequencies: Pi = — = 26.49 + 


0.07 days and P 2 = -b = 26.92 ± 0.07 days (lMassi & Jaronl 
120131) . The period l\ agrees with the value of 26.4960 ± 
0.0028 days ( Gregor y 12002b associated with the orbital period 
of the binary system and with the predicted periodical accretion 
peak. Period P 2 agrees well with the estimate by radio astrome¬ 
try of 27-28 days for the precession period. 

Recently, a Lomb-Scargle analysis of Fermi- LAT data 
around apastron revealed the sa me periodicities Pi = 26.48 + 
0.08 d and P 2 ^ = 26.99 + 0.08 d (ijaron & Massi 120141) . The sim¬ 
ilar behaviours of the emission at high (GeV) and low (radio) 
energy would imply that the emission is caused by the same 
population of electrons in a precessing jet. This important re¬ 
sult must be confirmed by simultaneous observations. GBI ob¬ 
servations were made from January 1994 until October 2000. 
Fermi-LAT has observed since August 2008. The fact that the 
two archives do not even overlap shows the persistence of vi and 
v 2 and characterizes them as permanent features of the system 
LS I +61°303. However, to investigate the electron population, 
one needs simultaneous monitorings. Lortunately, a new radio 
long-term monitoring is available. The Owens Valley Radio Ob¬ 
servatory (OVRO) 40 m telescope monitoring at 15 GHz began 


in March 2009. In this Letter we present a timing analysis of 
these data. Section 2 describes the data reduction and the results. 
Section 3 presents our conclusions. 


2. Data analysis and results 

2.1. OVRO observations 

Regular 15 GHz observations of LS I +61°303 from 54908.8 
MJD (March 2009) to 56795.0 MJD (May 2014) were carried 
o ut app roximately twice per week using the OVRO telescope 
(IRichards et al Jf20TTI) . The OVRO 40 m uses off-axis dual-beam 
optics and a cryogenic high electron mobility transistor (HEMT) 
low-noise amplifier with a 15.0 GHz centre frequency and 3 GHz 
bandwidth. The two sky beams are Dicke-switched using the 
off-source beam as a reference, and the source is alternated be¬ 
tween the two beams in an ON-ON fashion to remove atmo¬ 
spheric and ground contamination. Calibration is achieved us¬ 
ing a temperature-stable diode noise source to remove receiver 
gain drifts, and the flux density_scale is derived from observa¬ 
tions of 3C 286 assuming the lBaars et alJdl977l) value of 3.44 Jy 
at 15.0 GHz. The systematic uncertainty of about 5% in the 
flux density scale is not included in the error bars. Complete 
details of the reduction and calibration procedure are found in 
IRichards et al.1 (1201 It) . OVRO data vs time, long-term, and or¬ 
bital phase are shown in Lig. 1 

2.2. Timing analysis 

To search for possible periodicities, we used the Lomb-Scargle 
method, which is very effic ient for irregularly sampled data 
dLombl 119761 : IScarglel 19821) . We used the algorithms of the 
UK Starlink software package, PLRIOL0. The statistical signifi¬ 
cance of a period is calculated in PERIOD following the method 
of Lisher randomization as outlined in iLinnel Nemec & Nemeq 
( 1985 ). The advantage of using a Monte Carlo- or randomiza¬ 
tion test is that it is distribution-free and is not constrained by any 
specific noise models (Poisson, Gaussian, etc.). The fundamental 
assumption is that if there is no periodic signal in the time series 
data, then the measured values are independent of their observa¬ 
tion times and could have occurred in any other order. One thou¬ 
sand randomized time series are formed, and the periodograms 
are calculated. The proportion of permutations that give a peak 
power higher than that of the original time series would then 
provide an estimate of p, the probability that for a given fre¬ 
quency window there is no periodic component present in the 

1 http://www.starlink.rl.ac.uk/ 
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August 2008 - November 2014 
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Fig. 2. Lomb-Scargle periodograms of three independent long-term 
monitorings of LS I +61°303. Top: OVRO at 15 GHz. Centre: Ferrni- 
LAT in the energy range 100 MeV to 300 GeV using only data from 
the orbital phase interval <X> = 0.5 - 1.0 fsee I Jaron & Massi! l2014h . 
Bottom|_GBI at 2 GHz (empty circles) and 8 GHz (filled squares) (see 
iMassi & Jaron lf20T3 l). In all periodograms there are two periodicities 
with a false-alarm probability of between 0.00 and 0.01 (see Sect. 2). 
The ratio (R) between the intensity of the two spetral features at Pi and 
P 2 is different in the different periodograms and is R=1.9 at 2 GHz and 
R= 1.5 at 8 GHz for GBI data, R= 1.4 at 15 GHz and R=1.3 for GeV 
data. 

data with this period. A derived period is defined as significant 
for p < 1%. a nd a marginally significant one for 1% < p < 10% 
(iLinnel Nemec & Nemed] 19851) . Figure 2 shows the result for 
OVRO data. There are two periodicities Provro = 26.5 + 0.1 d 
and Pi, ovro = 26.9 + 0.1 d. The periods are both significant: they 
result in the randomization test with a false-alarm probability of 

p < 1%. 


2.3. Fermi-LAT 

The timing analysis in I Jaron & Massi ( 2014 ) for the whole data 
set gave Pi as significant (p < 1%) in the randomization tests 
even if it is a rather weak feature in the periodograms of their 
Fig. 3 a, b, and c. In other words, the timing analysis performed 
on the whole Fermi-data set yields that P 2 is significant even 
if the spectral feature is much lower than that at Pi. When the 
timing analysis is performed for emission in the orbital phase 
0=0.5-1.0, as described in iJaron & Massi 1 (120141) . the spectral 
feature at P 2 is comparable with that at Pi, but the randomization 
tests find a probability of false detection of p = 4%. The y- ray 
data used in that analysis did span the time period MJD 54683 
(August 05, 2008) to MJD 56838 (June 30, 2014). Now we have 
five more months of Fermi-LAT observations^'I’he data_ reduc¬ 
tion was performed as reported in IJaron & Massi (Imi . The 
periodogram for data in the orbital phase 0=0.5-1.0 shown in 
Fig. 2 has the same periodicities as in IJaron & Massi d2014l) , 
but with the important difference that now the randomization test 
gives both frequencies as significant, with a false-alarm proba¬ 
bility of p < 1%. 

2.4. 385 days of GBI data vs RATAN results 

The frequency resolution in a periodogram is related to the in¬ 
verse of the overall time interval of observations. The resolu¬ 
tion for GBI of Av = 0.0001 d 1 covers six times the difference 
in frequency between iq = 0.03775 ± 0.00010 d -1 and V 2 = 
0.03715 ± O.OOOlOd -1 determined bv IMassi & Jaron I (120 1 3l) . 
Indeed, the two spectral features are evident in the GBI peri- 
odogra m shown here in Fig. 2. Very recently, iTrushkin et a l.! 
( 2014 ) reported that the timing analysis of RATAN-600 radio 
telescope data did not find the two frequencies. The reported 
daily monitoring with the RATAN-600 radio telescope includes 
the time from 16 November 2013 (MJD 56612) to 6 Decem¬ 
ber 2014 (MJD 56997), that is, 385 days. This corresponds to 
a frequency resolution of 0.00065 d 1 that is even higher than 
the difference in frequency between v\ and V 2 of 0.00060 d 1 . 
We tested the insufficient frequency resolution by performing 
the timing analysis on 385 days of GBI data. The resulting pe¬ 
riodogram, sho wn in Fig. [3 is well comparable with that of 
ITrushkin et al. I (120141) . The two frequencies iq and V 2 are plot¬ 
ted in Fig. 3 by arrows. The low spectral resolution is still unable 
to resolve the frequency separation. 

3. Conclusions and discussion 

The spectral resolution in the timing analysis of OVRO data is 
able to distinguish the two spectral features found in GBI and 
Fermi-LAT data (Table 1). 

These two periodicities are physically related to the period¬ 
ical ejection (Pi) of electrons in a particular orbital phase to¬ 
ward apastron and to the precession (P 2 ) of the jet. The beating 
of these two periodicities probably causes the so-called long¬ 
term modulation. The peak flux density of the periodical radio 
outbu r st exhibits in fact a m odulation of 1667+8 d dGregorvI 
2002). Massi & Jaron I (1201 3l) have shown the consistency be¬ 
tween the long-term modulation and the beating of the two fre¬ 
quencies determined in the GBI timing analysis, thatis , Ph P at = 
P ~ 1 - P ~ 1 ) 1 (Table 1). Moreover. lMassi & Torricelfi-Ciamporfil 
20141) have shown that a physical model for LS I +61°303 of 
synchrotron emission from a precessing (P 2 ) jet, periodically 
(Pi) refilled with relativistic particles, produces a maximum 
when the jet electron density is at its maximum and the ap- 
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g*:? 

*£ i+ 
> 

Pi ± A Pi 
days 

Ratio 

^beat — A.Pbeat 

days 

Monitorings 

Duration 

years 

Reference 

26.5 + 0.1 

26.9 + 0.1 

1.4 

1782 + 630 

OVRO (15 GHz) 

5.2 

this work 

26.49 ± 0.07 

26.92 ± 0.07 

1.9 

1658 + 382 

GBI (2 GHz) 

6.7 

a 

26.49 ± 0.07 

26.92 ± 0.07 

1.5 

1658 + 382 

GBI (8 GHz) 

6.7 

a 

26.48 ± 0.08 

26.99 ± 0.08 

1.3 

1401 ±311 

Fermi-LAT (0.1 - 300 GeV), O = 0.5 - 1.0 

6.3 

b 


Table 1. Timing analysis results. The ratio between the intensity of the two spectral features at Pi and P 2 in the periodograms of Fig. 2 is given 
in the third column. Pbeat is define d as Pb eat = (P 7 1 - P, 1 ) -1 . Column six reports the duration of the related monitoring. The references are a: 
iMassi & Jaron I J201 3l> and b: Ijaron & Massi I d2014j) . 


proaching jet forms the smallest possible angle with the line of 
sight. This coincidence of the highest number of emitting parti¬ 
cles and the strongest Doppler boosting of their emission occurs 
with a period of (P~' - P 7 1 ) \ that is, Pbeat , or the long-term 
periodicity. 

The open question now is: why is the precession so close 
to the orbital period? While the micr oquasar SS433 and several 
known precessing X-ray binaries riLarwoodl 19981) have a preces¬ 
sion period longer by an order of magnitude than the orbital pe¬ 
riod (as predicted for a tidally forced precession by the compan¬ 
ion star), the different c ase of the microquasar G RO J1655-40 
was discovered in 1995 ( iHiellming & Rupenlll995i). This syst em 
has an orbital period of 2.601+0.027 days (iBailvn et al .11 19951) . 
IHiellming & Runenl d 19951) discovered a ra dio jet in this object 
with a precessional period of 3.0+0.2 days (Hiellming & Rupen 
H9951) . LS I +61°303 seems to be the second case of this new 
class of objects with orbital and precession periods that are rather 
close to each other. 


CD 

O 

Q_ 


March 1995-April 1996 



Fig. 3. Lomb-Scargle periodograms of 385 da ys of GBI radio data at 
2 GHz to be compared with the periodogram in lTrushkin et al. I d2014h 
also of 385 days. The spectral resolution of only 0.00065 d -1 cannot 
resolve the separation of 0.0006 d -1 between Vi = 0.03775±0.00010 d -1 
and v 2 = 0.03715 ± 0.00010 d -1 determined by IMassi & Jaron I ( 120131) 
with 6.7 yr of GBI data. 
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